Localized states due to expulsion of resonant impurity levels from the continuum in 

bilayer graphene 



3 



V. V. Mkhitaryan and E. G. Mishchenko 
Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA 

Anderson impurity problem is considered for a graphene bilayer subject to a gap-opening bias. In- 
gap localized states are produced even when the impurity level overlaps with the continuum of band 
electrons. The effect depends strongly on the polarity of the applied bias as long as hybridization 
with the impurity occurs within a single layer. For an impurity level inside the conduction band 
a positive bias creates the new localized in-gap state. A negative bias does not produce the same 
result and leads to a simple broadening of the impurity level. The implications for transport are 
discussed including a possibility of gate-controlled Kondo effect. 
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Introduction. Electronic properties of graphene are 
affected by short-range scattcrers (for the review see 
Ref. Remarkably, these properties significantly de- 
pend on the character of disorder Q ■ It has been demon- 
strated that the limit of strong scattering is realized via 
resonances occurring for vacancies [3|-|5( and adsorbed 
atoms Such resonances lead to the modification of 

the bare spectrum near the charge neutrality point and 
strongly affect the transport properties. 

Moreover, the localized impurity levels in a graphene 
host could potentially result in a formation of local mag- 
This situation is of a special in- 



ilHlil 



nctic moments 
terest as it opens up a perspective of controlled switching 
of the local moments by gating. As far as the transport 
properties are concerned the local magnetic moments fa- 
mously lead to the Kondo effect. The band structure of 
the monolayer graphene makes this effect different than 
that in normal metals 16-18|. On one hand, narrow- 
ing of the impurity level due to the vanishing density 
of states (DOS) in graphene, p{Ef) oc |-Ef| q , facilitates 
the formation of a local moment near the charge neu- 
trality point. On the other hand, the low DOS disfavors 
the Kondo resonance and results in the suppression of 
the Kondo temperature when Ep — > 0. In case of non- 
interacting electrons, a = I, the suppression is power- 
law. If the effects of electron-electron renormalization of 
the spectrum are included, so that a < I, the Kondo 
temperature is suppressed exponentially. 

The theory of resonant impurity states in bilayer 
graphene has received much less attention despite its 
greater technological potential. The most important 
property of bilayer graphene to this extent is its gap tun- 
able by external bias [19l42fj . It is therefore important to 



elucidate how strongly the gap is affected by the disorder, 
which generally tends to smear the gap. For short-range 
impurity potential there exist "tail-states" with DOS de- 
caying exponentially from the band edge [22 - 24 1 . Numer- 
ical studies also indicate the onset of well-defined midgap 
peaks in the presence of resonant impurities fioj ]. 

In the present paper we analyze the effect of the An- 





FIG. f : a) Bilayer graphene subject to the gap-opening bias 
2(7. The resonant impurity (black circle) is hybridized with 
7r-electrons of a carbon atom in the upper layer. Tunneling 71 
is indicated by the vertical dashed lines, b) Band structure of 
the system for different impurity energies eo- The left picture 
corresponds to a zero hybridization amplitude V . The right 
one illustrates the modification of the spectrum for finite cou- 
plings V and positive bias U > 0. The level with eo > U 
(blue) broadens due to coupling to the band electrons while 
generating an additional in-gap state. The level originally po- 
sitioned inside the gap (shown with green for a specific case 
eo = 0) remains localized but is renormalized towards lower 
energies. States overlapping with the valence band £0 < — U 
(red) are broadened without generating any additional in-gap 
states. The level at the top of the valence band (orange) leads 
to a square-root singularity in the density of states. 



derson impurity hybridized with the band electrons of bi- 
layer graphene host. We demonstrate that a single impu- 
rity level with energy £0 induces a localized in-gap state. 
This is a rather trivial statement whenever £0 fells within 
the (gate-induced) bandgap. More interestingly, we find 
that the in-gap state appears even if the impurity level 
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sits inside the continuum of itinerant states, i.e. above 
or below the bandgap, depending on the polarity of the 
applied interlayer bias. We consider a realistic bilayer 
setup, shown in Fig. 1, of an impurity sitting on top of 
the upper layer. Consider the case of a positive applied 
bias 2U/e as shown in the picture (— e is the electron 
charge). We emphasize that in this case the localized 
state appears no matter how high the impurity level Eq 
is above the charge neutrality point. To the contrary, 
the impurity state overlapping with the valence band, 
e < — U, docs not generate a localized state. It simply 
acquires a finite width (and the energy rcnormalization) 
from the interaction with the valence band states. 

The localized states of the similar nature that appear 
due to hybridization with an impurity level are known 
to occur in various systems. In a one-dimensional band 
with a finite bandwidth a pair of states split off near both 
band edges (2(| . In conventional s-wave superconductors 
similar bound states are long known 271 128|. However, 
they lie extremely close to the superconducting gap (typ- 
ically within ~ 10~ 3 of its width), and thus are hardly 
discernable against the divergent superconducting DOS. 
In contrast, in bilayer graphene it is possible to "place" 
the localized state virtually anywhere within the bandgap 
simply by controlling the applied bias. In particular, this 
can lead to the emergence of a local magnetic moment 
and a formation of the Kondo resonance, more robust 
than in a monolayer, even when the bare impurity level 
£o is above the Fermi level. 

Anderson impurity in bilayer graphene. The low- 
energy Hamiltonian of bilayer graphene is a 4 x 4 matrix 
in layer / sublattice space [29( . If in addition the energies 
of interest are less than the interlayer tunneling rate the 
Hamiltonian could be reduced further to the 2x2 ma- 
trix acting on the layer index. In the presence of a gap- 
opening electrostatic potential difference 2U such Hamil- 
tonian has the form [l[ [2^] , 



55 <« 



p 2 y ]<J x + 2p x p y a y ) +Ua z 



(1) 



where <7j are the Pauli matrices, and p is the electron 
quasimomentum. The effective mass m = 71 /2v 2 = 
0.03mo is determined by the interlayer coupling, 71 = 
0.35 eV, and the Dirac velocity, v = 10 6 m/s. This 
Hamiltonian acts on the two-component wavefunction, 



(i/'ipj il>2p)> w ffh VVp referring to the amplitudes 



on upper (a = 1) and lower (er = 2) layers. The to- 
tal Hamiltonian of the system in the second quantization 
representation is 

H = ^ P #(p)4 + e dU + Vj2(4pd + dtyip), (2) 



where the second term describes the localized level and 
the last term defines the hybridization of this level with 
band-electrons in the upper layer only. This choice, while 



not of principal importance for the calculations, ensures 
the asymmetry with respect to the polarity of U as we 
show below. The properties of the localized electrons are 
most conveniently described by its Green's functions, 



D(t) = -i(Td(t)d!(0)). 



(3) 



This function satisfies the equation which follows from 
Eq. ([2]). In the energy representation, 



{E-so)D(E) = l + Vg 1 {0,E). 



(4) 



The mixed Green's function in the latter equation is de- 
fined according to 

e ff (r,t) = -»Wa(r,t)d t (0)) (5) 

and satisfies the equation 

(E - H) aa ,Q a , (r, E) = VS(r)5 al D(E). (6) 

The two equations (j4|) and (J6j) form a closed system that 
readily yields the solution, 



D(E) 



E-e -V*G$(E) 



(7) 



where we introduced the Green's function of free band 
electrons, G^(E) = *£ p [E - H(p) + ir) signE]' 1 . In 
what follows we are mostly interested in the gapped re- 
gion, — U < E < U, where after simple calculations [3(J 
we obtain G ( ^(E) = -f Remarkably, in this 

region there exists a bound state which emerges as a sim- 
ple pole of the Green's function D(E), Eq. ([7]). Its energy 
satisfies the relation 



(8) 



Here we introduced the parameter, T = mV 2 /2, whose 
meaning is the inverse lifetime of the original impurity 
level high above the gap, E 3> U. Remarkably, the equa- 
tion ([8j is asymmetric with respect to the sign of the 
applied bias U. Let us concentrate on the case of posi- 
tive bias U > below. The generalization to the opposite 
case is almost evident and will be discussed later. It is 
easy to see that Eq. © has a true in-gap bound state 
solution, e, for any impurity level such that £0 > —U. 
Fig. [TJ illustrates the energy spectrum of the system for 
various values of sq. The most trivial case is realized 
when £0 falls within the bandgap — U < £0 < U. Here 
the energy of the resulting localized state is lowered by 
the hybridization, e < Eq. Moreover, impurity level in- 
side the positive energy band £0 > U always creates an 
in-gap localized state, regardless of how large £0 and/or 
how narrow the gap is. 

The effect of the impurity level is two-fold. First, its 
interaction with band electrons broadens its density of 
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FIG. 2: a) Dimensionless density of states (scaled with the 
coupling strength F) of the d-electrons plotted in units of 1 /V 
for r = 7i/2, £o = O.771, and three different values of bias 
U: zero bias (black dashed line), positive bias (red), negative 
bias (green). The bound in-gap state appears as a sharp peak 
for U > 0, but not for U < 0. Inset: development of a square- 
root singularity when the impurity level is approaching the 
top of the valence band from the inside of it. 



states. This could be readily seen from Eq. (0 with 
\E\ > U, where the self energy of the d-state becomes 



E+U 



The density 



purely imaginary, (E) - -. T vgLga . 
of states of the impurity level pd{E) is modified from the 
simple Lorcntzian, which is the case when U = 0: 



Pd{E) 



TVE 2 -U 2 /i 



(E~e ) 2 \E-U\+T 2 \E + U\ 



(9) 

Fig. [2] shows the density of impurity states for different 
values of the gap. The second feature is the appearance 
of the delta-function peak, corresponding to the new lo- 
calized state, Eq. ©, also indicated in Fig. [5] The total 
area under any curve is J pd{E)dE = 1. In evaluating 
this area one notices that the in-gap peak in Pd{E) has 
a spectral weight smaller than unity: 



Pd{E) 



(U -g) 3 / 2 (t/ + e) 1 ' 2 

{u -?) 3 / 2 (u + ?) 1 /2 + ur 



5{E-e). (10) 



The remaining spectral weight that completes the to- 
tal weight of the in-gap state to the unity comes from 



the band electrons. This can be verified by calculat- 
ing the band electron Green's function, G aa i (r, r', t) = 
—i(Tip a (r,t)ij}l,(r',0)), accomplished by constructing a 
system of coupled equations similar to Eqs. (|4]) and ((6]). 

Impurity levels overlapping with the valence band. 
When the impurity level resides inside the lower contin- 
uum of band states, £0 < —U, the broadening of the level 
is not accompanied by the creation of an in-gap state, see 
Fig. [TJ Far from the band edge — £0 U the lineshape 
is a simple Lorentzian but is distorted as £0 moves closer 
to the edge. At £0 = — U the density of states develops a 
singularity pd ~ T^ 1 \/U /\E + U\. Both these cases are 
shown in the inset on Fig. [2] 

The asymmetry between the positive and negative en- 
ergies £0 and the sensitivity to the sign of the external 
bias is a manifestation of breaking of the symmetry be- 
tween the two layers by the impurity position. If the 
impurity particle is hybridized equally with both mono- 
layers the symmetry is restored and the in-gap bound 
state is always present. It is found from an equation sim- 
ilar to Eq. ([5]), where the square root is replaced with 
2E/\/U 2 - E 2 . 

4x4 bilayer Hamiltonian. In the previous discussion 
we utilized the simplified two-band model, which is usu- 
ally justified when the typical energies are small com- 
pared with the interlayer tunneling, E, U <C 71. We are 
now going to discuss the new qualitative features that 
the inclusion of all four 7r-electron bands leads to. The 
full four-band Hamiltonian is a 4 x 4 matrix 1291 . 



ff(p) = 



/ U v P+ \ 

vp- U 71 

71 — U vp_ 

\ vp + -U ) 



[11) 



where p± — p x ± ip y . The first two rows/columns of Eq. 
(jlip refer to the two sublatticcs in the upper graphenc 
layer and the last two to the lower layer. The formal ex- 
pression for the total Hamiltonian of the system, Eq. ([lj , 
remains unchanged with the exception of the dimension- 
ality of spinors ip p . The effect of hybridization is the 
strongest when the impurity sits on a carbon atom be- 
longing to the first sublattice - the one not directly cou- 
pled by tunneling 7! to the second layer. A rather pecu- 
liar situation of an impurity residing on the other sublat- 
tice of the first layer is discussed below. 

The pole of Eq. ([7|) still defines the bound state with 
the corresponding self-energy given by the integral of the 
Green's function over momentum, 



^(Q)^_ U ~ E i 1 2 (U + E) + 2UE(U-E) 7T E 2 + U 2 
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One of the consequences of the expanded Halimtonian 
(ITTjl is that the true gap U = 7if//V7i + ^U 2 is found 
at a finite momentum. Since U and is slightly less than 
U, the logarithmic term in Eq. (Ti"2"]) is cut-off and repre- 
sents a non-significant renormalization of e . The most 
significant difference obtained in the 4x4 model is the 
behavior at E ~ — U, where Gii(E) changes sign and 
therefore allows the in-gap state solution expelled from 
the inside of the valence band, £o < — U, situation for- 
bidden by Eq. ©. Such states are residing very close to 
the top of the band as the overall effect is rather weak. 
For the binding energy we obtain the expression, 



U 



&r 2 u 5 



it(\s \-uy 



(13) 



applicable as long as |eo| — U ^> TU 2 /j 2 . 

Let us now turn to the impurity sitting on top of 
an atom belonging to the second (directly coupled by 
tunneling to the bottom layer) sublatticc. In this case 

the self-energy in Eq. ([7]) is given by (E) instead of 
G^(E). The former is determined by expression similar 
to Eq. (fT2")) where the term -f 2 (U + E) in the numerator 
of the second term is set to zero. There are two notable 
consequences of this change. First, all bound solutions 
for the impurity overlapping with the conduction band, 
£o > U, completely disappear. The only non-trivial states 
are those emerging from the valence band. Second, such 
states are expelled stronger than in the case considered in 
the previous paragraph. This is reflected in the increased 
binding energy, which is given by Eq. (fl3|) multiplied by 
the extra factor 4. 

Gate- controlled Kondo effect via formation of a local 
magnetic moment. A direct inclusion of the spin degree 
of freedom in the Hamiltonian © does not change the 
above physical picture unless the electron-electron inter- 
action is taken into account. In the presence of the spin 
degeneracy, this interaction is often described by the on- 
site Hubbard repulsion, which under some conditions can 
cause a magnetic instability [31j . Such a description, 
however, might not be adequate for the in-gap states 
found above. Indeed, as evident from Eq. (fTTJj) . a sig- 
nificant fraction of the spectral weight comes from the 
evanescent band states, in particular, when e is close to 
the band edges, ±U. We consider the most intriguing sit- 
uation of an impurity level high in the conduction band 
Eq ^s> U, so that the in-gap state is close to U . To esti- 
mate the magnitude of the Coulomb charging energy, we 
note that the radius of this state is R ~ [2m(U — e)] -1 / 2 . 
In turn, the binding energy can be found from Eq. ([8]) 
to be U - e ~ 2T 2 U/{e - U) 2 . Therefore the charging 



energy is E c 



iU/eq. The condition for 



the formation of a well-defined magnetic moment in the 
?i-doped bilayer (with the Fermi energy Ep > U) is that 
the double occupancy of the in-gap state is forbidden, 
E c > Ef — U. In terms of the gap-opening bias, this 



happens when U > Epe\j EbT 2 , where we introduce the 
Bohr energy, Eb ~ 1-47 cV. Under these conditions the 
gate-controlled Kondo effect can potentially be observed 
in bilayer graphene, where one has the additional benefit 
of a finite DOS, in contrast to a monolayer. 

Conclusion. An important difference between the 
above found resonant bound states with the previously 
predicted bound states due to short-range impurities 
[22l . 23 1 or vacancies [25| is that in our case the emerging 



states can reside anywhere inside the gap, not just close 
to the band edges. Wide range for the hybridization ener- 
gies for a number of impurities (H, CH3, OH, F) has been 
reported from the first principles from W ~ 0.7 eV in 
Ref. i to W ~ 5 eV in Ref. 0. Within the effective Dirac 
description this corresponds to V ~ (3 3 ^ 4 a/2h)W , where 
the interatomic distance a = 1.42 A. Using this estimate, 
T ~ 10 — 250 meV. Due to this large value of T, the depth 
of the induced in-gap bound state is comparable with U 
for values of £0 as large as a few V. In addition to the 
above discussed Kondo effect, the potential experimental 
implications include sub-gap gate-dependent optical ab- 
sorption, as well as the variable-range hopping transport 
via the impurity band, which would form in the presence 
of a finite density of impurities. 
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